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ABSTRACT 

The determination of velocities of stars from precise Doppler measurements 
is described here using relativistic theory of astronomical reference frames so as 
to determine the Keplerian and post-Keplerian parameters of binary systems. 
Seven reference frames are introduced: (i) proper frame of a particle emitting 
light, (ii) the star-centered reference frame, (iii) barycentric frame of the binary, 
(iv) barycentric frame of the Galaxy, (v) barycentric frame of the Solar system, 
(vi) geocentric frame, and (vii) topocentric frame of observer at the Earth. We 
apply successive Lorentz transformations and the relativistic equation of light 
propagation to establish the exact treatment of Doppler effect in binary systems 
both in special and general relativity theories. As a result, the Doppler shift is 
a sum of (1) linear in c -1 terms, which include the ordinary Doppler effect and 
its variation due to the secular radial acceleration of the binary with respect to 
observer; (2) terms proportional to c~ 2 , which include the contributions from the 
quadratic Doppler effect caused by the relative motion of binary star with respect 
to the Solar system, motion of the particle emitting light and diurnal rotational 
motion of observer, orbital motion of the star around the binary's barycenter, and 
orbital motion of the Earth; and (3) terms proportional to c -2 , which include the 
contributions from redshifts due to gravitational fields of the star, star's compan- 
ion, Galaxy, Solar system, and the Earth. After parameterization of the binary's 
orbit we find that the presence of periodically changing terms in the Doppler 
schift enables us disentangling different terms and measuring, along with the 
well known Keplerian parameters of the binary, four additional post-Keplerian 
parameters, which characterize: (i) the relativistic advance of the periastron; 
(ii) a combination of the quadratic Doppler and gravitational shifts associated 
with the orbital motion of the primary relative to the binary's barycentre and 
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with the companion's gravitational field, respectively; (iii) the amplitude of the 
'gravitational lensing' contribution to the Doppler shift; and (iv) the usual incli- 
nation angle of the binary's orbit, i . We briefly discuss feasibility of practical 
implementation of these theoretical results, which crucially depends on further 
progress in the technique of precision Doppler measurements. 

Subject headings: Gravitation - binaries: general - pulsars: general - interfer- 
ometry: optical 
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1. Introduction 

Binary stars represent perhaps the most valuable targets for stellar astrophysics. They 
have been a source of insight into the structure and evolution of stars, theory of radiative 
transfer, stellar magneto- and hydrodynamics, and the Newtonian theory of gravity, to men- 
tion just a few major topics (Shore 1994). After discovery of the first binary pulsar 1913+16, 
these objects have become excellent gravitational laboratories for testing general relativity 
by using radio observations (Taylor 1992, 1994). Nevertheless, optical observations of bi- 
nary stars continue to remain one of the most important sources of getting qualitatively new 
astronomical information. The reason is that the amount of binary stars observable at op- 
tical wavelenghts overwhelmingly exceeds the number of objects accessible for radio, X-ray, 
and/or 7-ray observations. Moreover, distribution of the relative orientation of binary orbits 
ranging from the edge-on to the face-on, allows one to obtain valuable information in study- 
ing different stellar phenomena. Therefore, increasing accuracy of optical measurements of 
binary stars is a real challenge for modern astronomy. This is why precision Doppler mea- 
surements of stellar spectra implemented recently for the search of extrasolar planets could 
open a new direction in binary star research. 

Traditional techniques in radial velocity measurements rarely achieve an accuracy better 
than about 200 m s -1 . Given such an uncertainty, usually the linear Doppler effect, i.e. the 
term of order of only v/c could be only measurable. Its measurement for the primary star of 
mass m s belonging to a binary system brings five classical Keplerian parameters of the stellar 
orbit. They are: projected semimajor axis a s , eccentricity e, orbital period P&, longitude of 
the periastron u, and the epoch T of the initial periastron passage. Combination of these 
parameters makes it possible to calculate the mass function for the binary system: 

m 3 c sm 3 i 

f{m s ,m c ) = - ■ -, 1 

[m s + m c ) 2 

where m c is the companion mass, and i is the inclination angle of the stellar orbit to the line 
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of sight. Obviously, the information which can be extracted from the linear Doppler shift 
alone is incomplete to determine all orbital parameters of the binary, including the masses 
of its constituent stars. Observations of additional relativistic effects are necessary. They 
include: relativistic advance of the periastron, quadratic Doppler and gravitational redshifts, 
monotonous decrease of orbital period due to the emission of gravitational waves, and effects 
of deflection and retardation of electromagnetic waves in the companion's gravitational field. 
All of them have been observed in binary pulsars (Taylor 1992) but only relativistic advance of 
periastron could be measured with the use of conventional spectroscopic technique (Semeniuk 
& Paczynski 1968, Guinan & Maloney 1985). 

Ozernoy (1997a, 1997b) pointed out that current accuracy of Precision Doppler Mea- 
surements (PDMs) using the iodine-based Doppler technique (Valenti, Butler & Marcy 1995, 
Cochran 1996) is able to catch the second-order in v/c effects. Moreover, he also has shown, 
by taking into account special relativity alone, that binary stars offer a unique opportunity to 
disentangle the linear and quadratic in v/c terms and extract such an important parameter 
as inclination angle. 

Even without addressing any concrete applications, a coherent, unambiguous interpre- 
tation of PDMs having an accuracy of ~1 m s _1 or better, requires an adequate development 
of a relativistic theory. By present, the basic principles to construct such a theory have been 
well established. Recently, they have been worked out in a series of publications by Brum- 
berg & Kopeikin (1989a, 1989b, 1990) (BK approach) and Damour, Soffel & Xu (1991, 1992, 
1993) (DSX approach). The main idea is to introduce one global and several local coordinate 
charts in a gravitating system consisting of N bodies to describe adequately the properties 
of space-time curvature both on the global scale and locally in the vicinity of each body. The 
subsequent application of the mathematical technique to match the asymptotic expansions 
of the metric tensor in different coordinate systems allows one to obtain general relativistic 
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transformations between these systems, which generalize Lorentz transformations of special 
relativity. Practical conclusions of two approaches are the same. In this paper, we use the 
BK approach to describe self-consistently the relativistic algorithm of PDMs. 

The paper is organized as follows. In section 2, seven appropriate coordinate frames 
are introduced. Sec. 3 deals with the Doppler effect in special relativity, which can be 
important for interpretation of PDMs in the situations, when effects of the gravitational field 
are negligibly small. The post-Newtonian coordinate transformations are sketched in section 
4, along with the derivation of the equation for propagation of photons in a gravitational 
field. The Doppler effect in general relativity is explored in Section 5. Parameterization 
of the Doppler effect and the explicit Doppler shift curve are given in section 6. Section 7 
outlines some observational implications of the theory. Finally, section 8 contains discussion 
and our conclusions. The approach developed in this paper was earlier reported in Kopeikin 
& Ozernoy (1996). 

2. Coordinate Frames 

A rigorous mathematical treatment of precision Doppler observations of a binary star 
requires the use of seven (4-dimensional) reference frames (RFs): 

(G)— barycentric reference frame of our Galaxy (cT,X) = (X°,X l ); 

(S 1 )— the Solar system's barycentric reference frame (ct,x) = (x°,x l ); 

(-B)— the binary system's barycentric reference frame (cs,z) = (z°,z l ); 

(C)— the star-centered reference framef] (cA,ff) = (77°, rj 1 ); 

5 It is important to emphasize here that A denotes coordinate time in the star-centered 
reference frame and is not a wavelenght of photon. 
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(E)— the Earth's (geocentric) reference frame (cu,w) = (w°,w' t ); 
(T)— topocentric reference frame of terrestrial observer (cr, £) = 

(P)— proper reference frame of a particle emitting light (cv,() — (C°)C*)-Here and 

hereinafter, the arrow above the letter denotes a spatial vector having three coordinates; 
index "0" relates to time coordinate, and small latin indices (such as i = 1,2,3) represent 
the spatial coordinates. 

The origin of each coordinate frame coincides with the center of mass (barycenter) of 
the respected system of gravitating bodies. For instance, the origin of the Solar system 
RF is at the center of mass of the Solar system, the origin of the emitting star's RF is at 
the center of mass of the star, and so on. The observer is regarded to be massless and 
placed at the origin of the topocentric RF. We assume that emission of light is produced by 
the atom placed at the origin of its own proper RF (P). Each RF has its own coordinate 
time. These times are related to each other by means of relativistic time transformations 
(Brumberg & Kopeikin 1989a, 1989b, 1990). It is worth noting that the coordinate time 
of the topocentric RF coincides precisely with the proper time of the observer measured by 
the atomic clocks, and the coordinate time of the RF of the atom emitting light coincides 
with its proper time. The barycentric RF of our Galaxy is considered to be asymptotically 
flat so that it covers all space-time. All other coordinate frames are not asymptotically flat, 
and they cover only restricted domains in space because of a non-zero space-time curvature. 
All coordinate systems are assumed to be nonrotating in the kinematical sense (Brumberg 
& Kopeikin 1989a). It means that spatial axes of all RF's are aligned and anchored to the 
outermost quasars whose proper motions are negligibly small. 

To derive equations describing the Doppler effect at the post-Newtonian level of accu- 
racy, we use the relativistic post-newtonian transformations between the coordinate frames. 
They have been formulated by Kopeikin (1988) and Brumberg & Kopeikin (1989a, 1989b), 
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and are discussed briefly in section 4. However, for pedagogical reasons, it is useful to con- 
sider the Doppler effect first in the framework of special relativity. The special relativistic 
approach is motivated by the fact that one can get the exact treatment of the problem under 
consideration, which will provide a guide of consistency for general relativistic calculations. 
It should be noted, however, that the special relativistic approach is not able to take into 
account the relativistic effects asssociated with an accelerated motion of bodies as well as in- 
fluence of gravitational field. Those efffects can be adequately considered only in framework 
of general relativity. 



3. Doppler effect in Special Relativity 

The Doppler effect in special relativity is usually considered only for two RFs: one is 
assumed to be at rest, and the other moves with respect to the first one with a constant 
velocity. Here we discuss a more realistic situation when five coordinate systems S,B,C,E 
and T, introduced in the previous section, need to be considered. This situation is rather close 
to the real astronomical practices and it might be applied to the interpretation of spectral 
observations of binary stars if the influence of gravitational fields could be neglected. In this 
section, we abandon, for the sake of simplicity, the reference frames G and P. The reason is 
that, for the moment, we want to avoid the discussion of terms caused by the motion of the 
Solar and binary systems about the center of our Galaxy, as well as motions of the emitting 
particles with respect to the star. Accounting for these effects will be done in later sections. 

The RF S is the basic one which is considered to be at rest. The origins of the RFs B and 
C are moving with respect to S with constant relative velocities Vb, and V^, respectively. 
The RF C is supposed to move with respect to the B with a constant relative velocity v c . 
Note that, in the relativistic approach, v c ^ *Vc- The RFs E and T move with constant 
velocities Ve and with respect to S. The relative velocity of the reference frame T with 
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respect to E is v T (again, v T ^ Y T ). 

Three different approaches for discussion of Doppler effect can be applied. They are 
based, accordingly, on the techniques of relativistic frequency transformation (Landau & 
Lifshitz 1951), successive Lorentz transformations (Weinberg 1972), and time transforma- 
tions along with the equation of light propagation (Tolman 1934, Brumberg 1972). 

3.1. Frequency transformation technique 

Let k a be the 4-vector of the electromagnetic wave propagating from the source of light 
to the observer. Here and hereinafter, greek indices run from to 3, i.e. k a = (k°, k l ). This 
is a null vector in the flat space-time. Therefore, k° = 2inj/c, and k l = —k°n l , where v is 
the frequency of the electromagnetic wave and the unit spatial vector n l is tangent to the 
trajectory of the light ray. For convenience, it is chosen to be directed from the observer 
toward the point of emission. Let u a = (u°, u % ) be the vector of 4- velocity of a massive 
particle. The time component is u° — I/7, and the spatial components are u l = u°/3 l , where 
7 = (1 — /3 2 ) 1//2 , is the (constant) Lorentz-factor, (3 l = v l /c, and v l is a spatial velocity of 
the particle. By contracting vectors k a and u a , one forms a scalar which is relativistically 
invariant as it is independent of the choice of reference frame: 

u a k a = u % k % - u°k° = invariant, (2) 

where the repeated spatial indices mean a summation from 1 to 3. 

Suppose that the source of light moves with respect to the observer with a constant 
speed v l . The 4- velocity of the observer in its proper RF is defined as u a = (1, 0, 0, 0), and 
the 4-velocity of the source of light is u a = 7 _1 (1,/?*). Let the frequency of the emitted 
electromagnetic wave be i>o, and the received frequency be v. Then, by applying equation 
(0) to the two different RF's one gets (boldface letters denote spatial vectors, and the dot in 
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between two spatial vectors stands for usual scalar product): 

z/o _ 1 + (ff-n) 

which is a well-known result (Landau & Lifshitz 1951) for the Doppler shift of the frequency 
of light emitted by a moving source and received by the observer at rest. Here the unit 
spatial vector n is measured with respect to the observer's RF. 

While applying this formula, we will use slightly different notations: v* for the frequency 
of the emitted light and v for the observed frequency, viz., u s , z/&,, and v e for the frequencies 
of light observed in reference frames S, B, and E, respectively. Let us also introduce a 
fractional frequency shift function, z = — — 1. Succesive application of equation (^) yields 

v*v h v B v e 1 + Q3 c -N fl ) 1 + Q9 B -N) (1-/31) 1 / 2 (l-ffi) 1 / 2 
v h v s v e v {l-m i/2 (1-/3I) 1/2 l + (/3 £ -N) l + ({3 T -N E y U 

Here N = (x* — x)/ |x* — x| is the spatial unit vector tangent to the light ray and having 
components measured with respect to the coordinate system S, and N#= (z* — z)/|z* — z| 
and = (w* — w)/^* — w| represent the same tangent vector with the components mea- 
sured relative to the systems B and E, respectively. Equation @ contains the dimensionless 
particle velocities (3 B = Vg/c, P c = vc/c, /3 E = Vf;/c, and /3 T = Vt/c. Coordinates with the 
asterisk concern the point of emission of light, and coordinates without asterisk are related 
to the point of observation. It is worth noting that the components of vectors N# and 
do not coincide with those of vector N because of the relativistic aberration of light. 

A remarkable feature of the formula @ is that it represents the Doppler effect as a 
product of four different multipliers. Each factor describes transformation of frequency of 
light from one reference frame to another. It would be straightforward to generalize this 
result for the description of the Doppler effect in the event of as many reference frames 
as necessary. Because of importance of equation (|), it is instructive to derive it by using 
different techniques and then to compare the results. 



- 11 - 



3.2. Lorentz transformation technique 

Lorentz transformations from the reference frame S to B is described by the matrix A^j 
with the components (Weinberg 1972, Brumberg 1972, 1991): 

Ag> = 7b 1 , Ag = Ag = -1b X Pb, 

(5) 

A l i = ^ + (7b 1 - IWPbPb- 
Similarly, the Lorentz transformation from the reference frame B to C is described by the 
matrix A£f with the components: 

A° c ° = 7c 1 , Ag = A§ = -7c^c, 

(6) 

Ag = + (7c 1 - l)Pd*Pc& 
The Lorentz transformations from the reference frame S to E, and from E to T are given 
respectively by the matrices A^f and A^with the components: 

A°° = 7b 1 , Ag = = -7^/3! , 

(7) 

Ag = ^ + (7b 1 - , 



(8) 

A" = ^ + (7^ - 1)/3t 2 ^4 • 

The relationship between time components of a light ray's 4- vector is given by the successive 
Lorentz transformations (the repeated greek indices imply summation from to 3): 

k° c = k^k^kl , (9) 
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k° T = h^h^kl , (10) 

where kg, k^, and kj> are components of the light ray vector referred to the RF's S, C, and 
T, respectively. 

By substituting the matrices of the Lorentz transformations into equations (PD,(PT|) and 
defining the null vector kg = 2ttvc~ 1 (1, — N l ), we obtain after straightforward calculations 
the following result: 

1+z ~ {l -^ )l ' 2 (1 "^ )1/2 x (11) 

1 + (ggjg + 7b(/3 c -N) + (gggg) + (1 - 7 B)^ 2 (/3 c -/3 j? )(/3 B -N) 
1 + 03 B .N) + 7s (/3 r " N) + (/3 T -/3 B ) + (1 - 7e)P¥ 2 (P t -Pe)(Pe-K) ' 

At first sight, it looks quite different compared to equation (f|). However, by making rela- 
tivistic transformation of vectors Ng and to vector N in equation (§), one can readily 
show that both expressions are completely identical. In the rest of thid section, we derive, 
for the reader's convenience, the transformation law between vectors and N. (The trans- 
formation law between vectors and N is obtained similarly by replacing index E for B 
and coordinates w l for z % .) 

The transformation between spatial coordinates of RFs S and E is given by (Weinberg 
1972, Brumberg 1972, 1991): 

w* = A%(x> - V 3 E t) , (12) 

where the transformation matrix is defined in equation (0). In its explicit form, the 
transformation fll2"[) reads 

1 



w = x — + 



- 1 



VWx-V^ 



V 2 
1 E 



Let us express the coordinates of radius- vector connecting points of emission and obser- 
vation in the RF E through the coordinates of vector N. From eq. ( [TBI ) and equation (BTJ) 
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for light propagation one gets: 



wi — w 



and, as a consequence, 



w* — w 



|x*-x|A*|(AP+/%) 



1 + (Pe ■ N) 



(14) 



1 (1 - V|/c 2 )V2 • 

Now it is easy to obtain the relationship between vectors and N, whisch is given by: 



(15) 



m 



(1-/3|) 1/2 
1 + 09js-N) 



N +1e Pe + \1e ~ l ) m 



01 



(16) 



Inversely, vector N l is obtained from equation ([16]) after replacements N l 
N\ and (5% — > 



AT 2 



: ! \i/2 



1 - 09*-N B ) 



£1 



(17) 



The transformations ([[!]) and ( |T7| ) represent, in fact, general expressions for the rela- 
tivistic aberration of light rays. This can be seen from the relativistic law of addition of 
velocities (Weinberg 1972, Brumberg 1972, 1991). In case under consideration, it is given 
by: 



(1-/3|) 1/2 
l + c-H/Vv) 



v +1e Pe + \1e - l ) & 

Pe 



where v l and V 1 are the relative velocities of a particle with respect to RFs E and S, 
respectively. For the light particle (photon) these velocities are v % = —cN E , and V % = —cN l . 
Having substituted them to eq. dl8|) , one obviously gets eq. (JjT 



Finally, using equations (]T6|) - (0), one obtains: 

l + (N E .f3 T ) 1 + (N.(3 E ) ' (19) 

One can see from equation ( |T9"D and a similar expression for 1 + (Nb-{3 c ) that equations (|j) 
and ( |TT| ) are identical. The advantage of eq. ( |TT]) over (Q) is that only one vector N enters 
eq. (|TT|), instead of three vectors N#, Ng, and N in eq. (0). 
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3.3. Time transformation technique 



Previous techniques used to derive the Doppler equation have not taken into account 
an essential fact of separation of the two events - emission and observation of light - in 
space-time. In fact, we have implicitly assumed that the null vector k a is the same at the 
points of emission and observation of light. However, this is only true for a very special case 
of negligent gravitational field and propagation of light in vacuum. In general, this condi- 
tions are not met. Therefore, a more advanced technique is required to tackle the Doppler 
effect appropriately. Such a technique, based on the integration of the equation for light 
propagation from the point of emission to the point of observation, establishes a relationship 
between coordinates of the 4- vector of a photon at these two events. Transformation laws of 
time scales between different RFs are to be taken into account as well. This approach, being 
rather general and straightforward, can be applied to analyse any particular situation. In 
this section, we consider time transformation technique in special relativity only. Its appli- 
cation to observations of binary stars in the framework of general relativity will be discussed 
in later sections. 

The equation of light propagation, in the absence of gravitational field and interstellar 
medium, is quite simple: 



where t* is the instant of photon emission and t is the instant of observation of the photon, 
both measured as coordinate time of RF S, in which x* = x(t*) is the point of emission, and 
x = x(i) is the point of observation. It is worth noting that, although the instants t and t* 
belong to the same RF S, their increments At and At* are different because of a relative 
motion of the source of light and the observer. 

One can see from equation (p0|) that when the influences of gravitational field and the 




(20) 



1 

c 



— X 



(21) 
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medium are both negligent, the components of the vector k a are constant everywhere on 
the light ray's trajectory. This makes clear why we do not care about a point in space-time 
in which earlier we calculated the Doppler shift. However, in a more general situation, the 
equations fl20|) , fl21 ) are not so simple, and this point has to be appropriately taken into 



consideration. 

Doppler effect is described by the function 1+z = ^f, where = 1/AA* is the frequency 

of the emitted light, and v = 1/ At is the observed frequency. By taking the time intervals 

to be infmitesimally small, we get a differential formula: 

dr du dt dt* ds* 
du dt dt* ds* d\* ' 



which is nothing more but a simple rule for differentiation of an hierarchical function /(A*) = 
-^(^(^(^(A*))))). This result demonstrates as well that the Doppler effect can be presented 
as a product of several multiplyers. A difference between equations ( |22| ) and (f|) is that in 
(f22|) we use coordinate times of the respected RFs and distinguish explicitly the points of 
emission and observation of light. Meanwhile in equation (^) only proper frequencies of the 
electromagnetic wave are considered. The advantage of eq. (^) is that, for its derivation, one 
needs to know only relativistic transformations between time scales, whereas transformation 
law between spatial coordinates is not required. As we shall see later on, this advantage is 
very helpful while tackling the Doppler effect in general relativity. 

To calculate the time derivatives at the points of emission and observation, one needs 
incorporating time components of the Lorentz transformations (H) - (§) between different 
RFs. They are: 

_ u-c-\(3 T -w) 

T ~ (l-/3f.)l/2 ' [Z6) 

_ U-c 1 {f3 B ■ x*) 

*~ (1-/32)1/2 > W 
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A - = * a -fr'-' • < 26 > 

Let us remind that observer is fixed with respect to the RF T, and the source of light is 
fixed with respect to C. Therefore, (3t = c _1 v^ = c~ l dtw /du, and (3 C = c~ l \c — c^dz^/du. 
Velocities of the observer and the source of light relative to the RF S are = dx/dt and 
Vc = dx.*/dt, respectively (it is important to note that Vj- 7^ Vr and Vc 7^ vc). Therefore 
one obtains from (^J) - ( |2"6T ) : 

^ = (1-^) V2 , (27) 
du 



du 1-((3 E - V r ) 



<K (1-/?1)V2 
dt* (1-/5I) 1/2 



(28) 
(29) 



In addition, differentiation of equation ( |2"T| ) gives: 

dt l + c-^N-Vc 



ds* l-C-^g-Vc) ' 

" S 1 - £) l2 - (30) 



(31) 



dt* l + c-!(N- V T ) ' 
Substitution of expressions fl2~7|) - (|31|) into eq. (^) gives for the Doppler shift: 

1 1 _ (1~/%) 1/2 (1-^) V2 1-c-^-Vt) l + c-HN-Vg) 

(1 - /?S)i/2 (l _ £2)1/2 ! _ c - 1(y g s . Vo ) 1 + c-HN • V T ) ■ 1 J 

This equation does not coincide apparently neither with (^), nor with ([IT]). Nethertheless, 

taking into account relativistic transformations between vectors Ng, Ny, and N as well as 

the law of addition of spatial velocities one can readily show that all three expressions for 

the Doppler effect are identical. Indeed, with the use of equation ([18]) it follows that 

l-c-H/yvr) _ (i-/?i) 1/2 

(1-/31)V2 1 + ((3 E -(3 T ) ' 

and 

1 + ((3 E -N) + T s(/VN) + ((3 T ■ (3 E ) + (1 - IeW^t ■ Pe)(Pe-N) 



(33) 



1 + cT 1 (N- V T ) 



1 + (Pe ■ Ar) 

(34) 
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These relationships, along with ones obtained from (|33| ) - ([34]) after replacement of indices 
T to C and E to B, allow to see that expression (|32|) for the Doppler effect coincides with 
(|TT1) and, consequently, with (|). 

This completes the derivation of the exact equations for the Doppler effect in special 
relativity. These equations could be expanded into powers of 1/c to get an approximate 
solution. Unfortunately, we would not be able to apply directly those expressions to real as- 
tronomical practices since gravitational fields of the Solar system, the binary system, and the 
Galaxy give contributions comparable with the special relativistic quadratic Doppler shift. 
Thus, in order to explore the Doppler effect in general relativity, it is important to elaborate 
approximative analytical methods . To tackle this problem, we apply the relativistic theory 
of astronomical reference frames developed by Kopeikin (1988) and Brumberg & Kopeikin 
(1989a, 1989b). 

4. Coordinate transformations in General Relativity 

Transformation laws between reference frames in general relativity generalize the Lorentz 
transformations of special relativity. They can be derived in two steps. First of all, the 
explicit form of metric tensor in different RFs are obtained by solving the Einstein equations 
with relevant boundary conditions. Then, the general relativistic transformations between 
the RFs are derived using the method of matched asymptotic technique. A clear and simple 
introduction to this technique is given in Brumberg & Kopeikin (1990). Here we give the 
transformation laws in the form which is suitable for discussion of the Doppler effect with a 
more than sufficient accuracy. To derive the Doppler shift, we apply the technique based on 
time transformations (Sec. 3.3). Thus, there is no need for development of relativistic part of 
space-time transformation between spatial coordinates, which will be given hereinafter only 
in the Newtonian approximation. 
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4.1. Transformation between topocentric and geocentric reference frames 

This transformation law is given by: 



1 

T = U 



+ 0(c" 4 ), (35) 



C = w i -w l T + 0(c- 2 ), (36) 

where u;^(«) and v^(u) = dw\jdu are geocentric spatial coordinates (RF E) and velocity 
of the observer, respectively; and $t is the geopotential at the observer's location site. It 
is worth noting that the quantity f|i>f + ^t) is constant on the geoid surface. In eq. (|35|), 
the tidal gravitational potential of external bodies is not included since it is negligibly small. 
Geocentric coordinates and observer's velocity both depend on time. Once the observer 
(spectrograph) is at the surface of the Earth, its and are precisely calculated using 
the data of the International Earth Rotation Service (IERS). If the observer is on board of 
a satellite, its motion can be derived using the satellite monitoring service. 



4.2. Transformation between geocentric and solar barycentric reference frames 



This transformation is found in the form: 

1 



1 

t - -r 



c 2 L 



/ (^l + ^e) dt + v k E {x k - x|)J + 0(c" 4 ) , (37) 

iu* = x* - 4 + 0(c' 2 ) , (38) 

where x|(t) and v^(t) = dx\jdt are respectively the spatial coordinates (RF S) and velocity 
of the geocentre relative to the barycenter of the Solar system; and Ue is the gravitational 
potential of the Solar system at the geocenter. If the external (with respect to the Earth) 
bodies of the Solar system are approximated by massive point particles, then 

= E ^ (39) 



k=l 
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where m k is mass of the body k; r k is the distance from the body k to the geocentre, and 
the sum is taken over all the external bodies of the Solar system. The potential <£> T is not 
included in Ue, in accordance with general principles of construction of relativistic theory 
of astronomical reference frames (Kopeikin 1988, Brumberg & Kopeikin 1989a, Brumberg 
& Kopeikin 1989b). The tidal gravitational potentials of the bodies external with respect 
to the Solar system are not included either, because they are too small to be important in 
the calculations of the Doppler effect. Barycentric coordinates and velocities of the Earth 
and other bodies of the Solar system can be calculated using the contemporary numerical 
theories of their motions (Standish 1982, 1993). 

4.3. Transformation between the solar and galactic reference frames 

This transformation law reads: 



where Xg(T) and V s fc (T) = dX^/dT are spatial coordinates and velocity of the barycentre of 
the Solar system with respect to the barycentre of our Galaxy; W$ is the gravitational po- 
tential of the Galaxy at the barycentre of the Solar system (the potentials $t and Ue should 
not be included). The galactic coordinates, velocity of the Solar system, and gravitational 
potential of the Galaxy at the Solar system barycentre are all not well known quantities 
so far. To measure them more accurately would be one of many practical implications of 
precision Doppler measurements of stars. 




(40) 



x i = X l -X l s + 0(cT 2 ), 



(41) 
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4.4. Transformation between the binary and galactic reference frames 



This transformation is similar to eq. ( pEUD and is given by: 




(42) 



x* = X i -X B + 0(c- 2 ), 



(43) 



where Xg(T) and V£(T) = dX^/dT are spatial coordinates and velocity of the barycentre of 
the binary system relative to the barycentre of our Galaxy, respectively; Wb is the gravita- 
tional potential of the Galaxy at the barycentre of the binary system (gravitational potential 
of the binary system should not be included). 

4.5. Transformation between the stellar and binary reference frames 

This transformation is similar to eq. ([57]) and has the form: 



where z c (s) and v c (T) = dz c /ds are spatial coordinates and velocity of the primary star 
relative to the barycentre of the binary system; and Uc is the gravitational potential of the 
companion star. Gravitational potential of the primary star should not be included in Uc 
for the same reason why the potential U-& does not include geopotential $x- The potential 
Uc is given by: 



where m c is the companion mass, and r is the distance between the two stars in the binary. 




(44) 



i7W-afc + 0(c- a ), 



(45) 



4.6. Transformation between the proper frame of an emitting atom and 

stellar reference frame 



This transformation law is given by: 




1 r 



/ 



it£ + $p) d\ + v*(r) k -4)\+0(c- 4 ), 

C = v l -v l P + o(c~ 2 ), 



(47) 



(48) 



where ?7p(A) and Vp(X) = drj^/dX are spatial coordinates and velocity of an emitting atom 
relative to the star, respectively; $p is gravitational potential of the star at the point of the 
atom's location. Obviously, coordinates and velocity of a single emitting atom cannot be 
determined since the integral flux of the stellar radiation is only observed. Motion of the 
atom and gravitational potential of the star both causs the broadening of spectral lines in the 
stellar spectrum. This unfortunately complicates the precise measurement of the Doppler 
shift. In order to simplify discussion of this problem as much as possible, we assume here 
that Vp, an average thermal velocity of atoms, is constant in time, and $p, gravitational 
potential of the star at the altitude of the spectral line formation, is also a constant. 

4.7. Time transformation between instants of emission and observation 

Time transformation between instants of light emission and observation is obtained from 
the solution of equation for light propagation in vacuum, which is described by the equation 
of isotropic geodesic line (Weinberg 1972). Solution of this equation has a simple form in the 
galactic reference frame G so that the time interval between the instants of light emission, 
T*, and observation, T (T > T*), is given by: 



where X\ are the galactic coordinates of the emitting atom at the instant of emission, and X % 
are the galactic coordinates of the observer at the instant of light observation. Relativistic 




X*-X|+A 5 (T,T*), 



(49) 
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correction As is of order of 0(c 3 ). It describes the Shapiro time delay (Shapiro 1964) in 
the gravitational field: 

^ C 2 J Gv(X a ^ |X^ — X a | "I - |X — X a I |X* — X| (50) 

a c 3 |X„, — X a | + |X — X a | — |X* — X| 
where subscript a stands for a body a that deflects light rays, X l a are its spatial coordinates 
taken at the moment T a of the closest approach of the photon to the body (Klioner & 
Kopeikin 1992, Kopeikin et al. 1998). It can be shown (Brumberg 1972) that the main 
term in the Shapiro delay depends logarithmically upon d, the impact parameter of the 
light ray (for more detail see also the paper (Kopeikin 1997)). The contribution to the 
Doppler shift caused by the Shapiro delay is proportional to (v/c)(r g /d), where v is the 
characteristic relative velocity, and r g = 2GM/ c 2 is the gravitational radius of the deflector. 
This estimate makes it obvious that the contribution floTf ) to the Doppler shift produced by 
the Shapiro delay can be only substantial in the nearly edge-on binary systems containing 
invisible relativistic companion - a neutron star or a black hole. Functions X(T) and X^T*) 
can be decomposed into a sum of vectors 

X(T) = X S (T) + x E (T) + w T (T) + 0(c" 2 ), (51) 

X*(T*) = X B (T,) + z c (T,) + n P (T m ) + 0(c- 2 ), (52) 

where the relativistic terms come from the relativistic part of the transformation of spatial 
coordinates. In subsequent calculations of the Doppler shift in general relativity, we will use 
eq. (|49|) coupled with these expansions. 



5. Doppler effect in General Relativity 

5.1. General equation 

Frequency of the emitted light is related to the proper time of the emitting atom as 
= 1/A-u*, where Av^ is the period of the emitted electromagnetic wave. Frequency of 
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the observed light is v = 1/Ar, where At is the period of the received electromagnetic 
wave. The Doppler shift z in frequency is calculated as a product of the appropriate time 
derivatives: 

^ dr dr du dt dT dT* ds* dA* 

di;* (it tff dT* (is* d\* dv* ' 

where dr/du, du/dt, dt/dT are taken at the instant of observation; dT*/ds*, ds*/d\*, d\*/dv* 
are taken at the instant of emission; and dT j dT* is calculated by finding the differential of the 



left and right hand sides of equation (|9|) for propagation of light. Thus, equation (|53|) is not 
just the usual time derivative taken at the same point of space-time. On the contrary, this 
is a two-point function that relates two events separated in space and time and connected 
by an isotropic worldline. 



5.2. Expansion into a series in 1/c 

Time derivatives at the point of observation are obtained by direct differentiation of 
equations 



(37), and (40), which describe relativistic transformations between different 
time scales in the Solar system. All these derivatives are taken at the point of observation: 

1 



dr _ l 1 

du c 2 



-v%< + $ T (w T ) 



0(c 



du 1 
~dt = ~ 
dt 1 



77^1 + ^e(xe) 



dT 



-V S 2 + W S (X S ) 



1 • 



4) + 0(c" 4 ), 



+ -Vg(X«-X$)+0(c 



(54) 

(55) 
(56) 



Here a| = dv^/dt is the acceleration of the geocentre relative to the barycentre of the Solar 
system, and Vj? = dV^/dT is the acceleration of the barycentre of the Solar system with 
respect to the barycentre of our Galaxy. 



Calculations of time derivatives at the point of emission can be done with the help of 
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equations ©, d§, and (0): 



i+ M^ +$p(r/p) ] +o(c " 4) ' 



(57) 



ds* 



1 + ? [2^ + ^ c(Zc) 



1 



-a^(^-4) + 0(c- 4 ), 



(58) 



c 



1 + ~ 2 \\v£ + W B (X B )1 - - 2 V^ - X k B ) + 0(c" 4 ). 



(59) 



(is* 



Here Oq = dvc/dt is the acceleration of the emitting star with respect to the binary system 



with respect to the barycentre of our Galaxy (all quantities are calculated at the point of 
emission). 

The function dT^/dT depends on two instants of time, viz., emission and observation 
of light. We have found that it is more convenient to transform dT^/dT to the instant of 
emission alone so that to express the final result through the instantaneous relative velocity 
of the barycentre of the binary with respect to the barycentre of the Solar system. It enables 
us to exclude from the final equation for the Doppler shift the poorly known velocities of the 
binary and the Solar systems with respect to the centre of mass of our Galaxy. To complete 
this procedure, we introduce the notations as follows: 

• R l = x B (t*) = X B (T*) — Xg(T*) + 0(c~ 2 ) is the relative distance between the Solar 
system and binary barycentres taken at the instant of emission; 

• K 1 = R l /R is the unit vector directed toward to the barycentre of the binary (this 
vector slowly changes due to proper motion fi; 

• v 1 = dR l /dt* is the relative velocity of the binary's barycentre relative to the barycentre 
of the Solar system, taken at the moment of emission; 

• v r = (K ' v)K % is the radial velocity of the binary's barycentre; 



barycentre, and V£ = dV^ 



/dT* is the acceleration of the barycentre of the binary system 
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tiy = [Kx [v x K]]* = fiR is the transverse velocity of the binary's barycentre. 



The two-point time derivative dT^/dT can be found by means of calculation of differen- 
tial of equations fli^p, (|5TD, and fl52"). This results in: 

dT l + c-^N- V*) + er 3 F* 



dT» l + c-^N- V) + c- 3 F ' 
(N-V*)-(n*-V*) + (n*-V a ) (N • V*)+(n* • V*) - (n* • VJ 



-R*a + Ra — D 



R*a + Ra + D 



F = 2Gj2M a 



(N-V)+(n-V)-(n-V_J + (N • V)-(n • V) + (n ■ V e 



(60) 
(61) 

(62) 



-R*a + R a — D R^ a + R a + D 

where V, V*, and V a are galactic velocities of the observer, source of light, and deflecting 
body a, respectively; R* a = |X*(T*) - X a (T a )|; R a = |X(T) - X a (T a )|; D = |X(T) - X a (T a )|; 



N 



n 



defined 


as: 


X* (X 1 * 


) - X(T) 


|X(T) 


-X a (T a )| 


X*(T*) 


- X a (T a ) 


IX* (T + ) 


-X a (T a )| 


X(T)- 


-X a (T a ) 


|X(T) - 


-X a (T a )| ' 



(63) 

(64) 
(65) 



Furthermore, equation (|60| ) is expanded into the powers of c l , and it can be simplified using 
the relationships: 



R* a + Ra — D 



d 2 R* a + Ra 



+ 0(d 4 



2 R* a FL a 
R* a + R a + D = 2(R W + R a ) + 0{d 2 ) 



n* = N 



R* 



» = " N + ;H N (£) + 



(66) 
(67) 
(68) 

(69) 



where £ = [N x [R* a x N]] = — [N x [R a x N]] is the vector of impact parameter d pointing 
from the deflector to the light ray: d = |£| -C min(i? a , i?* a ). 
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If we only take into account the Shapiro effect in the binary system, then R a 3> 
and equation (|60"D takes the form: 



^ = l + -(N-V,)-i(N-V) + + i(N-V) 2 -- 2 

dl* c c <r c z 



1 + - (N • V*) - - (N • V) + + — (N • V) - — (N ■ V*) (N • V) + 

(70) 



W -4 + ^(v,-v c) . 



We have neglected in ( fT0|) all terms of the order of V 3 /c 3 and higher, as well as those 
'mixed' terms from the differentiation of As, which are of the order of (R# a / R a )(r g /d)(V/c), 
(dj 'R* a ) 2 {r g l 'd)(V '/ c), and so on, where V is the characteristic relative velocity of the com- 
panion relative to the primary, r g = 2GM C / c 2 is gravitational radius of companion, and d is 
the impact parameter of the light ray. 

To continue, we expand the function Xg(T) in eq. fl5T| ) into the time series near the 
instant T* : 

X(T) = X s (i;) + V S (T»)(T - T*) + x E (T) + w T (T) + O [(T - T,) 2 ] , (71) 

and, instead of T — T*, we substitute the r.h.s. of equation (f49"D- The result is used to expand 
the unit vector N from eq. (^) into the powers of parallactic terms of the order of Xe/R, 
Zq/R, and so on. One gets: 

N = K + 7T C - 7r E -i [Kx [V S (T„) x K]] + 0( e - 2 ) + O^V) + 0(tt 2 ). (72) 

c 

Here the term depending on the velocity Vs describes the secular aberration. The binary 
orbital parallax tzq (caused by the orbital motion of the star) as well as the annual parallax 
7Te (caused by the orbital motion of the Earth) are given by: 

7r c =i[Kx[z c xK]], (73) 
7r E =i[Kx[x E xK]]. (74) 
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In equation ([72]) ) we also neglect, as currently unmeasurable, all terms of the order oiw^/R 
and rjp/R. Similarly to decomposition of functions X(T) and X*(T#), given by eqs. ([51]) and 
(|52|), decomposition of velocities V(T) and V*(T*) can be done. Then it is straightforward 
to show that 



V*(T*) - V(T) = v + v P (T») + v c (T*) - v E (T) - v T (T)— V s (T,)i? + 0(c" 2 ). (75) 

c 



5.3. The Doppler shift 

After substitution of the intermediate equations of the previous subsection into the basic 
equation (|53|), the final result for the Doppler shift takes the form: 



z(t) = z c + z RQ + z EQ - z R - z e - z s - z M , 



(76) 



where partial contributions are given by: 



zc = 1 c vr + ? (|4 + \^ 2 R 2 + W s - W B - (K • V S )R) + 
+± (\vl + $ P ) - £ (|4 + $ T ) , 



(77) 



z m = — (K ■ v T ) - -(1 + c-%)(K • v E ) 

c c 



(78) 



2* + Lfe - (K • v^ 3 



pi 

- constant - 



+ U E - (K • v E ) 



. periodic 



(79) 



z R = — (K • v P ) - -(1 + c~V)(K • v c ), 
c c 



(80) 
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v f + u c ) -My + f/o ) ... • (81 > 

^ / constant u V ^ / periodic 

^ = ^^ln[r*-(K •!•«)], (82) 

z M = ^(K.vc)(K.v E )-p(/i-Uc)- ( 83 ) 

Here fx — k — V T /R is the vector of proper motion of the binary's barycentre; r R is the 
radius- vector of the primary star relative to its companion; = |r^|; and we neglect all 
terms of the order of 10~ 10 and higher. The nature of the partial components in eq. (76) is 
as follows: 

The term zc contains a linear Doppler shift caused by the radial velocity of the emitting 
particle. It also includes both the quadratic Doppler and the gravitational shifts caused by 
the relative motion of the binary and the gravitational potential of our Galaxy, respectively. 
Contribution from (|i>t + ^t) causes broadening spectral lines in the primary's spectrum. 
The geopotential term (j^vp + $p) is constant in time, and all temporal variations of Zc are 
expected to be caused by a radial acceleration of the binary and/or its proper motion. 

The term zr q describes a linear Doppler shift caused by the rotational motion of the 
terrestrial observer with velocity v T and the orbital motion of the Earth's centre of mass 
with velocity vg. This term includes the radial component of the binary's relative velocity. 

The term ze® includes a sum of quadratic Doppler and gravitational shifts caused, 
respectively, by the orbital motion of the Earth relative to the barycentre of the Solar system 
and the gravitational fields of the Sun and planets. 

The term zr describes a linear Doppler shift caused by the radial velocity of source 
of light vp relative to the star's centre and the radial component of the orbital velocity of 
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the star's centre of mass Vq- This term, like zrq, includes the relative radial velocity of the 
binary. 

The term ze is a sum of quadratic Doppler and gravitational shifts caused by the 
orbital motion of the primary star with respect to the barycentre of the binary and by the 
gravitational field of the companion. 

The term zs represents a Doppler shift caused by the Shapiro delay in propagation of 
light in the companion's gravitational field. This effect can only be detected in the nearly 
edge-on bynary systems. Its magnitude is generally negligibly small. 

The term zm describes a Doppler shift caused by the effect of coupling of motions of 
the Earth and the primary star. 

6. The explicit Doppler shift curve 

6.1. The necessity of parameterization 

Equation fl7B| ) as such cannot be used for reduction of observational data. It should be 
re-written in a way which would clearly pinpoint the measurable parameters. We have also 
to assign the proper instant of time ("exposure mid-time" t) to any particular observation 
of stellar spectrum (Cochran 1996). Moreover, the exposure mid-time should be properly 
referred to the instant of light emission. The determination of the exposure mid-time is a 
rather difficult technical problem and we do not discuss it here (see Cochran 1996). As for 
the relationship between the exposure mid-time and the instant of light emission, it follows 
from the equation of light propagation P9] ) and has the well-known form extensively used, 
e.g., in pulsar timing data reduction programs (Taylor & Weisberg 1989, Doroshenko & 
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Kopeikin 1995): 

(1 + zc)K = t-t + A RG + A EQ + A R + A E + A s . (84) 

Here to is the initial epoch of observations; A^ Q , A EQ , A R , A E , As are respectively the Romer 
and Einstein delays in the Solar system; and A R , A E , A s are, accordingly, the Romer, Ein- 
stein, and Shapiro delays in the binary system. Their explicit expressions can be found in 
Damour & Taylor (1992), Taylor & Weisberg (1989), and Doroshenko & Kopeikin (1990, 
1995). 



6.2. Convenient vectors for tracking the binary system 

Let us introduce a triad of the unit vectors (Io, Jo, K) attached to the barycentre of the 
binary system (see Fig.l). The vector K is directed from the Solar system barycentre toward 
that of the binary system, and vectors Io, Jo are m the plane of the sky with I directed to 
the east, and J to the north celestial pole. Two other sets of the unit vectors, (I, J, K) and 
(i, j, k) are also introduced, which are related to (I , Jo, K) by means of two spatial rotations 
(Damour & Deruelle, 1986b): 

I = cos fl I + sin fl J , i = I , 

J = — sin I + cos Jo , j = cos i J + sin % K , (85) 

K = K , k = — sin % J + cos % K . 

In the above transformations, the angles fl (0 <fl < 2n) and i (0 < i < it) designate the 
longitude of the ascending node of the primary's orbit and the inclination of the orbit to 
the plane of the sky, respectively. Vector I is directed to the ascending node of the binary's 
orbit, and vectors (i, j) lie in the orbital plane in the sense of orbital motion. 

Vector K slowly changes due to a proper motion of the binary 

K = K + /i(t-* ), (86) 
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where t is the current time, and to is the initial epoch of observations. Therefore, the relative 
velocity of the binary's barycentre with respect to the Solar system is given by: 

V = R{fi a l + fi s J ) + v R K + v R (t - t ), (87) 

where R is the distance between the binary and the Solar systems; Vr is the relative radial 
velocity (vr = R) at the initial epoch to] vr is the radial acceleration; \i a and /i<5 are the 
respective components of the proper motion of the star in the sky. 



6.3. Relativistic terms in the Doppler shift 

Relativistic perturbations of the orbit of a binary system are described in Klioner & 
Kopeikin (1994). Using the D amour- Deruelle relativistic parameterization of the orbital 
motion (Damour & Deruelle 1985, see also Klioner & Kopeikin 1994), we get with the 
necessary accuracy: 

-(Ko-vc) = K s [cos(u + A) + ecosu} + 0(c" 3 ), (88) 
c 

K s = nx s (l-e 2 )- l/2 , (89) 

where x s = a s sin i/c is the projection of the semimajor axis a s of the primary's orbit onto 
the line of sight; and n = 2ir/Pb is angular frequency of the orbital motion (P& being the 
orbital period) given by 



n 



( 



GM 



V a R 



1/2 



m p m c „\ GM 



1+ 



(90) 



M 2 J 2a R c 2 . 

Here m s and m c are masses of the primary star and its companion, respectively, M = m s +m c , 
o>r — cisi^s + rric)/ m c + 0{c~ 2 ) is the semimajor axis of the primary's relative orbit in 
harmonic coordinates (Damour & Deruelle 1986b), and e is the eccentricity of this orbit. 
The angle u in eq. fl88D is the longitude of periastron, which includes a contribution of its 
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relativistic advance: 

lo = lu + kA e , (91) 

where u>o is the position of the periastron at the initial epoch, and k is the post-Keplerian 
parameter of relativistic advance of the periastron (Robertson 1938, Damour & Schafer 1985, 
Kopeikin & Potapov 1994): 



c 2 a R (l - e 2 ' 



As an example, for a binary system with the parameters m s = 3M , m c = 1.4M , and 
e = 0.4, the magnitude of the relativistic advance is about 0°.l yr~ l , 1.2" yr _1 , and 0.004" 
yr^ 1 for the orbits which semimajor axes are 10 12 cm, 10 13 cm, and 10 14 cm, respectively. 

The angle A e entering eq. ( pi]) is the eccentric anomaly, which is related to the time 
through the true anomaly U and the third Kepler's law: 

l + e\ 1 ' 2 u 



A P = 2 arctan 



tan — 



(93) 



.1-eJ 2 
U - e sin U = nX + a, (94) 

where a is the (constant) orbital phase at the epoch of the first passage of the periastron. 
Additional Newtonian perturbations of the binary orbit (whenever they are observationally 
important) may be included into equation (|88|) using the usual approach based on the orbital 
osculating elements (e.g. Shore 1992, p. 34). 

Coupling of orbital and proper motions of the binary gives the term 
1 K s 

-(fj. ■ v c ) = — \{n a cos Q + fis sinQ)S(u) + cosz(/x Q sin Q — fis cos Q)C(U)] , (95) 

c sin i 

where Q is the longitude of the ascending node of the orbit, and functions C(u) and S(u) 
are given by 

C(U) = cos(cj + Ae) + e cos u , (96) 
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S(U) = sin(^ + A e ) + e sin 00 = _^ZM. (97) 

Cud 

The quadratic Doppler effect plus gravitational shift of the frequency in the companion's 
gravitational field are given by 

Vn t t 1 Gm c \2m s + m c (l — e 2 )] Gm c \m s + 2m c ] e . , , 

-f + Uc = -r , 2V \ ^ f-z 5 cos A e . 98 

2 2 a fl (l - e z )(m s + m c ) a R (m s + m c ) \-e l 

Comparision of equations ( |9"8|) and (|8"ID yields 

Gm^ T 

(^constant = 2c2(Xr ^ + m j + 7 > (") 

(ze) periodic = TcosA e , (100) 

where a new relativistic post-Keplerian parameter T is given by 

Gm c \m s + 2mJ e 

T= 2 ) ° v i a- 101 

c z aR[m s + m c J 1 — e z 

For the parameters of binary systems given below eq. ©, the magnitude of T is about 
1.3 • 10~ 7 , 1.3 • 10" 8 , and 1.3 • 10~ 9 , respectively. 

Finally, for the "gravitational lens" term zs we obtain 

S{esinv4 e — sini [cos(o; + A e ) + ecoscj]} 

Zs = 1 : — 7 TTTi 2vT72 : T> ( 102 J 

1 — ecosu — sins [sin u; (cos w — e) + (1 — e z ) L/z coscusmwj 

where the third relativistic post-Keplerian parameter S is defined as 

_ 2Gm c n 

*~ C 3 (l_ e 2)l/2" ^ 

It is worth noting that the inclination angle % defines the shape of the function ( |102j ). There- 
fore, it can be considered as the forth post-Keplerian parameter (Taylor 1992). 

The effect of gravitational lensing is, under usual circumstances, rather small and it will 
be problematic to measure it. For instance, if the binary system consists of a main sequence 
star m s = 3M & and a relativistic companion m c = 1.4M Q on the orbit having the relative 
semi-major axis cir = 10 12 cm and eccentricity e = 0, the magnitude of zs is only 1.1 • 10~ 9 , 
2.4- 10 -9 , and 7.7 ■ 10 -9 for the orbital inclinations sin 2 = 0.95, 0.99, and 0.999, respectively. 
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7. Implications of the Doppler shift curve 

The set of equations derived in the previous section makes it possible to analyse PDM 
observations of binary stars on a quantitative basis. In this section, we explore how to 
disentangle various effects entering the basic equation (76) for the Doppler shift and extract 
measurable parameters. 

7.1. Effects of Earth rotation and orbital motion 

In order to extract a pure effect caused by the primary's motions, the Earth rotation and 
orbital motion have to be subtracted form the total Doppler shift. It can be easily done using 
machine readable data on the Earth rotation parameters (IERS Annual Report), the Earth 
spatial coordinates and the geocentre's velocity (Standish 1982, 1993) as well as position 
and proper motion of the binary star taken from the astrometric catalogue. Since this paper 
only deals with principal topics, we do not develop here an exact technical framework for 
calculating zrq, zeq and consider them in the following as well predictable functions of time. 

7.2. Effects of constant part of gravitational field and relative motion between 

the Solar system and the binary star 

The effect of relative motion between the Solar system and the binary star, has the main 
contribution of the order of 0(c _1 ), is associated with the radial velocity of the binary's 
barycentre and its radial acceleration. Terms of the order of 0(c -2 ) include the transversal 
velocity component fj,R squared and the acceleration of the Solar system's barycentre relative 
to the barycentre of the Galaxy V s . In addition, the constant parts of the gravitational fields 
of our Galaxy, the Solar system, the binary system, and the geopotential as well as the 
quadratic Doppler shift caused by the orbital motion of the Earth and the primary star, all 
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contribute at the level of 0(c~ 2 ). The geopotential, the quadratic Doppler shift caused by 
the orbital motion of the Earth, and the constant part of the gravitational field of the Solar 
system on the Earth's orbit all can be calculated based on the gravimetric data and modern 
ephemerides with an accuracy which high enough to exclude those terms from the function 
z. The remaining terms can be used to extract an additional information on the distribution 
of gravitational field in our Galaxy. 

7.3. Effects of the orbital motion of the primary star and proper motion of 

the binary system 

The classical Doppler effect associated with the orbital motion of a binary system is 
well known. It allows to measure five Keplerian parameters: (i) n, the frequency of the 
orbital motion; (ii) a, the initial orbital phase; (iii) e, eccentricity; (iv) uu , initial position of 
the periastron; and (v) x s = a s smi/c, the projected semimajor axis of the primary's orbit. 
Among these five parameters, two parameters, n and x, while combined together, make it 
possible to determine the mass function of the binary, f(m s ,m c ). In the binary with an 
invisible (compact) companion, the knowledge of mass function sets up an upper limit to 
the mass of the companion. In the event when the Doppler shift curve of the companion is 
also observable, this limit can be put even tighter. 

Measurements of relativistic effects in the orbital motion can provide a unique tool to 
determine the masses of stars in the binary without bias. In binary pulsars, such a procedure 
is widely used to determine the masses of neutron stars (Taylor 1992). 

The proper motion of the binary leads to a gradual secular change in the observable or- 
bital elements x s and u. This situation is quite similar to that in the pulsar timing (Kopeikin 
1994, 1996; Arzoumanian et al. 1996). Indeed, one can see from equations (|95|) - fl97|) that, 
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because of the smallness of proper motion, the term -(/j, ■ v c )(t — t ) can be entirely absorbed 
into -(Kq-Vq) by means of redefinition of parameters x s and oo. As a result, the observable 
values x° bs and uo ohs are shifted from their physically meaningful values x s and oo by 

x f s = Xs + 6x s , oo obs = 00 + 600, (104) 

where 

5x s = x s cot i {—^a sin Q + fis cos fl) (t — t Q ), (105) 

and 

5oo = esc i (fi a cosQ + /i^sinfi) (t — to)- (106) 

It is important to emphasize that the parameter x s changes because of secular variation 
of the inclination angle i due to proper motion. Meanwhile the semimajor axis a s remains 
constant, because proper motion does not cause any dynamical force acting on the orbital 
plane. Hence, equation (|105|) can be re- written in a form similar to ( |106|) : 

5i = (—fi a sin Q + fis cos Q) (t — t ). (107) 

The increments 5x s and 5oo depend on time linearly and can appear in observations as 
small secular variations of the Keplerian parameters x p and oo. Specifically, one gets: 

5x s = 1.54 • 10 _16 x s cot i (—Ha sinfl + fi$ cosfi) [s s ], (108) 

Soo = 2.78 • 10~ 7 esc i (jj a cos + n$ sin f2) [deg yr^ 1 ], (109) 

where the values x s ,fi a , and fi$ are expressed in seconds (s) and milliarcseconds per year 
(mas/yr), respectively. 

It is worth noting that, for binary systems with a negligibly small eccentricity, the effect 
of a secular variation in oo (along with the relativistic advance of periastron) is absorbed by 
the re-definition of the orbital frequency and therefore, in such systems, it is not observable 
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at all. Indeed, whenever the eccentricity is negligibly small, the argument uj + A e takes the 
form 



2tt 

DC 

b 

lobs 



LU + A e = u + —^(t-t ), (110) 



where P b is the observable value of the orbital period: 

(111) 



pobs = p b 



Pb 

1 — k — — esc i(fj, a cos Q + fig sin Q) 

Z7T 



Here Pb is the physical value of the orbital period and we neglected in ( |1 1 1|) all terms 
nonlinear in small parameters. 

If both the classical and relativistic perturbations of orbital motion are negligibly small, 
then as can be seen from eqs. ( 1Q8|) and ( [109| ), the observable secular variations of x s and 



uj parameters could be used to determine both the ascending node of the binary's orbit and 
the inclination angle of the orbit, i. 

Finally, it is worth noting that the term "§(/it • v c ) entering the function zm has the 
structure similar to the term -(ft ■ v c )(t — to). This does not include an explicit dependence 
on time and, therefore, only leads to a constant shift of the orbital parameters x p and uj, 
which therefore cannot be determined. 



7.4. The post-Keplerian parameters: k, T, and sini 

The post-Keplerian parameters k,T,$s, and sini can be measured in binary systems 
having relativistic orbits. Of the set of these parameters, the parameter T contributes at 
the highest order 0(c~ 2 ). However, it can be only disentangled with difficulty from the 
classical parameter K s cosu, similarly to what happens in binary pulsars (Brumberg et al. 
1975). The separation is possible only if the relativistic advance of the periastron, k, is high 
enough to measure the change in K s cosu. It should be noted that the term i^ecosu; that 



-38- 



enters eq. (88) is constant in the systems with negligibly small orbital perturbations. It has 
a secular change as the parameter oj is not constant. 

The parameter k contributes only at the level of 0(c~ 3 ) but in a secular way. This makes 
it measurement both easy and accurate, which has been done in a number of photometric 
and spectroscopic binaries (see, e.g., Shakura 1985, Khaliullin 1985). If parameters T and 
k are both measurable, then along with the mass function this would allow to determine 
separately the masses of both stars and the orbital inclination. If the classical perturbation 
of parameter u (for instance, caused by the oblateness of stars) is also substantial, then 
observations of parameter T would allow to separate the relativistic contribution to the 
advance of periastron from the classical one. This could be used to infer the oblatenesses of 
the stars. 

The two remaining parameters S and sini contribute, in a quasi-periodic way, at the 
level of 0(c~ 3 ). If it is done independently of the other parameters, S and sin 2 can only 
be measured in the nearly edge-on binary systems via the determination of the amplitude 
and the shape of function zs, in a manner similar to the pulsar timing (Taylor 1992). The 
range parameter S defines the amplitude of zc, and sini characterizes its shape. Once the 
parameters $5 and sini are determined, the masses of both stars can be obtained with the 
use of the mass function. 

8. Conclusions and Discussion 

For the reader's convenience, we summarize here the main conclusions of this paper, 
along with the references to relevant equations. 

1. As discussed in Introduction, Precision Doppler Measurements (PDMs), which pro- 
vide accuracy better than a few meters per second, measure more than just the radial com- 
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ponent of the velocity - they also catch a contribution of the transverse component, i.e. the 
terms of the second order in v/c. A source with perioically changing velocity components, 
such as a binary star, allows a disentangling different velocities and extracting additional 
(post-keplerian) parameters of the binary. To this end, a detailed relativistic theory of the 
Doppler shift is required. 

2. In special relativity, the Doppler shift of a spectral frequency from a binary is given 
by equation (4) or by equivalent expressions (11) and (32). The calculated Doppler shift 
includes the contributions from: (i) the motion of the binary's barycentre in the Galaxy; (ii) 
the motion of the primary star in the binary; (iii) the motion of the Solar system barycentre 
in the Galaxy; and (iv) the Earth motion relative to the Solar system's barycentre. 

3. In general relativity, accounting for additional effects is necessary, which includes, 
among others, gravitational field in the binary and acceleration of the primary star relative 
to the binary's barycentre. The total Doppler shift is given by equation (76), which partial 
components are presented by equations (77)-(83), with detailed comments about the physical 
meaning of each term given at the end of Sec. 5. 

4. Presence of periodically changing terms in equation (76) enables us disentangling 
different terms and measuring, along with the well known Keplerian parameters of the bi- 
nary, four additional post-Keplerian parameters (k, T, and smi) as well. The first three 
of them are given by equations (92), (101), and (103), respectively. The k parameter charac- 
terizes the relativistic advance of the periastron; T characterizes the quadratic Doppler and 
gravitational shifts associated with the orbital motion of the primary relative to the binary's 
barycentre and with the companion's gravitational field, respectively; 3 characterizes the 
amplitude of the 'gravitational lensing' contribution to the Doppler shift given by equation 
(102); and % is the usual inclination angle of the binary's orbit (the value of % defines the 
shape of the 'gravitational lensing' term). 
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5. The post-Keplerian parameters k,T,5s, and sin 2 can be measured in binary sys- 
tems, which are sufficiently close so as to make the relativistic effects measurable (on the 
other hand, as discussed below, for too close binaries with very short periods and therefore 
with very high orbital velocities there is a potential problem with the determination of the 
exposure mid-time while performing PDMs). 

Feasibility of practical implementation of the theory developed in this paper, crucially 
depends on further progress in PDM techniques. Analyzing the spectrum of a star with 
respect to the 'velocity metric' based on the spectrum of the iodine absorption cell allows 
one to remove most of the wavelength drifts of the spectrograph and the detector. It can 
provide an ultimate precision in measuring radial velocity of a star attaining 1 m s _1 or 
better (Cochran 1996). It is interesting to compare this precision with that in millisecond 
pulsars timing technique. Available data indicate that the limiting accuracy of pulsar timing 
measurements, 5t, with present techniques is a few microseconds, or less, over timespan of 
many years (Taylor 1992). Uncertainty in the velocity measurement is the product of radial 
acceleration of the body under consideration and the error in timing measurement. The star's 
radial acceleration with respect to the binary's barycentre is proportional to 4n 2 a sin i / P b 2 ~ 
2nv/Pb, where v is the radial component of orbital velocity, a is the semimajor axis of the 
star's orbit, and i is inclination angle of the orbit. Thus, for the typical binary pulsar PSR 
B1913+16, where 5t = 15 /is, Pb ~ 28000 s, and the ratio v/c ~10~ 3 , one gets the precision of 
timing velocity measurements Sv ~ 2irc(v/c)(5t/P b ) ~ 0.01 m s _1 . Thus, the PDMs of binary 
stars based on iodine absorption cell technique cannot currently be considered competitive 
with timing technique for binary pulsars as concerned the precise tests of relativistic gravity. 
Nevertheless, relative accuracy of PDMs (Sv/c ~ 3 • 10~ 9 ) is comparable with the magnitude 
of the second or (in some cases) even third order for relativistic perturbations in binary 
systems. Therefore, implementation of relativistic theory for proper tackling with such 
precision measurements is inevitable. 
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The prospects for PDMs combined with the relativistic theory of Doppler shift presented 
above seem to be even more bright if one takes into consideration that a new generation of 
instruments for PDMs is currently emerging. Among them is the use of Fourier Transform 
Spectroscopy with the Navy Prototype Optical Spectrometer being built by the U.S. Naval 
Observatory (Armstrong et al. 1998), which is projected to achieve in a near future a velocity 
resolution of only 0.3 m/s (Hajian 1998). 

An accuracy with which the orbital parameters of a binary can be determined by apply- 
ing the PDMs is restricted by an uncertainty related to the determination of the exposure 
mid-time. It is interesting to compare PDMs with pulsar timing where one can actually 
measure the arrival time of signals and not only the dopier shift. This allows a phase- 
connected solution for astrometric, spin, and orbital parameters of the pulsar which contains 
a fit to integer numbers. That is the reason for the high accuracy in pulsar timing exper- 
iment. More simply, the precision, with which the position of a pulsar on its orbit can be 
determined, is given by a relationship (Sr) PT ~ c • (St) PT . Let us assume, for convenience, 
that PDM gives an infinite precision in determination of stellar velocity but the exposure 
mid-time is determined with an error (St) PDM , which is about 1 s (Cochran 1996). Then 
inaccuracy in determination of the star's orbital position is given by (Sr) PDM ~ v ■ (St) PDM 
~ 2na sin i/P b (St) PDM . A comparision of the two expressions above gives: 

(5r) PT P b (St) PT ' 1 ] 

where x = asini/c. For binary systems having orbital parameters like those in PSR 1913+16 
with 2nx/Pb ~ v/c ~ 1CT 3 and (St) PT ~ 15 /is one has (Sr) PDM = 67 (Sr) PT ~ 300 km. 
However, it is worth emphasizing that since the accuracy of timing measurements has a low 
limit about 1 fis one can achieve a better determination of the orbits of binary systems in 
which v/c < 10~ 6 , i.e. for the systems with very long orbital periods. In this sense, PDMs 
are much better suited to search for planets orbiting the extrasolar stars than pulsar timing. 
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In this paper, we have only considered the post-Newtonian theory for PDMs of binaries 
consisting of a primary (optical) star and an (invisible) compact companion. Although in 
many binaries, especially spectroscopic ones, measuring the Doppler effect for both stars 
would provide additional possibilities for disentangling the orbital parameters, this kind 
of binary seems to be less appropriate for PDM measurements. Indeed, in close binaries, 
where relativistic effects could be measurable with PDMs, those effects would be severely 
contaminated by tidal interaction of the components and stellar winds, whereas in wide 
binaries the post-Newtonian terms are expected to be rather weak. 

It is worthwhile to remind that normal stars subjected to the PDMs, in most cases, 
cannot be considered as point masses, unlike neutron stars or black holes. Due to this 
reason, the classical perturbations will presumably be the most important sources of orbital 
parameters' variations (Shore 1994). However, even in the situation when the classical 
perturbations dominate, measuring of (or proper taking into account) the relativistic effects 
will serve as a tool to better understand the nature of the process(es) responsible for the 
orbital parameters' variations. In this respect, interesting observational targets for PDMs 
might be massive main-sequence stars with radio pulsars in binary systems like PSR B1259- 
63 (Wex et al. 1988) or PSR J0045-7319. Timing observations reveal (Lai, Bildstein & 
Kaspi 1995) that it is possible to measure orbital evolution caused by hydrodynamical effects 
associated with the optical companion but induced by the tidal gravitational field of the 
companion (the pulsar) and/or intrinsic rotational motion of the primary star. PDMs of 
such a system, to be done complementary to timing observations, would allow to measure 
tidal oscillations of the optical star and therefore to essentially improve our knowledge about 
these systems. Another interesting application of PDMs could be an anomalous binary 
systems like DI Her, where some discrepancy was found between the prediction based on 
general relativity and the observed motion of the periastron (Guinan & Maloney 1985). It 
also can be explained by the dynamical influence of tidally-generated oscillations on the 
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orbital motion of stars but accuracy of the employed observational technique was not good 
enough to test this hypothesis. 
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Fig. 1. — Angles and orientation conventions relating the orbit of the binary system to the 
observer's coordinate system and the line of sight. The orbital plane is inclined at angle % 
with respect to the plane of the sky. The angle Q is the longitude of the ascending node of 
the orbital plane. 



